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We report on a matched-filter search for gravitational wave bursts from cosmic string cusps using 
LIGO data from the fourth science run (S4) which took place in February and March 2005. No 
gravitational waves were detected in 14.9 days of data from times when all three LIGO detectors 
were operating. We interpret the result in terms of a frequentist upper limit on the rate of gravita- 
tional wave bursts and use the limits on the rate to constrain the parameter space (string tension, 
reconnection probability, and loop sizes) of cosmic string models. 
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I. INTRODUCTION 

Cosmic strings are one-dimensional topological defects 
that can form during phase transitions in the early uni- 
verse jU [2]. Topological defect formation is generic in 
grand unified theories (GUTs), and cosmic string produc- 
tion specifically is generic in supersymmetric GUTs [3J. 
In string theory motivated cosmological models, cosmic 
strings may also form (and are referred to as cosmic su- 
perstrings to differentiate them from strings formed in 
phase transitions) |H O E [TJ El [H EHl EH E21 [JS] . Cosmic 
strings and superstrings may produce a variety of as- 
trophysical signatures including gamma ray bursts |15j . 
ultra-high energy cosmic rays [2], magnetogenesis [25] . 
microlensing, strong and weak lensing [TBI H71 HH1 fHfl 120"] . 
radio bursts [21], effects on the cosmic 21 cm power spec- 
trum [22], effects on the cosmic microwave background 
(CMB) at small angular scales [2U [25], and effects on 
the CMB polarization [2J]. 

Cosmic strings and superstrings can also produce pow- 
erful bursts of gravitational waves [27] [25] |29"1 |3U1 I3T]. 
The most potent bursts are produced at regions of string 
called cusps which acquire large Lorentz boosts. The 
formation of cusps on cosmic strings is generic and cusp 
gravitational waveforms are simple and robust [531 ES]. 
The large mass per unit length of cosmic strings com- 
bined with the large Lorentz boost may result in signals 
detectable by Earth-based interferometric gravitational 
wave detectors such as LIGO [32] and Virgo [33] . Thus, 



gravitational waves may provide a powerful probe of early 
universe physics. 

The LIGO detector network is comprised of three laser 
interferometers. Two of them are located at the Hanford, 
WA site: a four-kilometer arm instrument referred to as 
HI, and a two- kilometer arm instrument referred to as 
H2. A second four-kilometer interferometer located at 
the Livingston, LA site, is referred to as LI. LIGO's 
fourth science run (S4) took place between February 22, 
2005 and March 23, 2005. The configuration of the LIGO 
instruments during the fourth science run (S4), is de- 
scribed in [36 . The sensitivity of this run was signifi- 
cantly better than that of previous runs: at the low fre- 
quencies relevant to this search, close to a factor of ten 
more sensitive than the previous science run S3, though 
still about a factor of two less sensitive than LIGO's most 
recent science run (S5), which was at design sensitivity. 

In this work we report on the results of a matched-filter 
search for bursts from cosmic string cusps performed on 
14.9 days of S4 data. We implement the data analysis 
methods described in |31j using a simple triple coinci- 
dence scheme. No gravitational waves were detected and 
we interpret the result in terms of a frequentist upper 
limit on the rate using the loudest event technique |37j . 
We use the upper limit on the rate to constrain the pa- 
rameter space of cosmic strings models. The sensitivity 
of the LIGO instruments during the S4 run does not al- 
low us to place constraints as tight as the indirect bounds 
from Big Bang Nucleosynthesis [3H]. In the future, how- 
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ever, we expect our sensitivity to surpass these limits for 
large areas of cosmic string model parameter space. 

In Sec.|ll]we discuss data selection, data analysis tech- 
niques, and describe the analysis pipeline. In Sec. |III| 
we describe the computation of the rate of accidental 
events we expect to survive the thresholds and consis- 
tency checks of the pipeline (the so-called background), 
and compare it to the events that made it to the end 
of the pipeline in our search (the so-called foreground). 
In Sec. |IV| we show how we estimate the sensitivity of 
the analysis using simulated gravitational-wave signals. 
We compute the efficiency of our pipeline, the fraction 
of simulated signals that we detect, as a function of the 
strength of the signals. In Sec. [V] we show how to es- 
timate the rate of burst events we expect, the effective 
rate, using the efficiency curves and the cosmological rate 
of events. We show the constraints our data place on the 
parameter space of cosmic string models. We conclude 
in Sec. ED 

II. DATA ANALYSIS 
A. Data selection and conditioning 

All available S4 science data when all three instruments 
were operating (triple-coincident data) were used except 
for periods 

1. with overflows in the error signal digitizer, 

2. when airplanes flew over the detector sites, 

3. thirty seconds prior to loss of lock (loss of reso- 
nance of the Fabry-Perot cavities in the arms) of 
any instrument, 

4. of excessive wind, 

5. of excessive seismic activity, and 

6. with calibration uncertainties larger than 10%. 

The total time of triple coincident data available after 
these cuts is 14.9 days. Calibration of the data used 
in this analysis was performed in the time domain |39j . 
The data were high-pass filtered near 30 Hz to remove 
un-necessary low frequency content, and down-sampled 
from the original LIGO sampling rate of 16384 Hz to 
4096 Hz. 



B. Matched-filters and templates 

For each of the three LIGO instruments we then per- 
formed a matched-filter search on this data for gravita- 
tional bursts from cosmic string cusps, i.e. linearly po- 
larized signals of the form [25] 

h+(f) = Br i/3 e(h - /) e(/ - /,). (i) 



The amplitude of the cusp waveform is B ~ 
GfiL 2 / 3 /(c 3 r), where G is Newton's constant, p is the 
mass per unit length of the string, L is the size of the 
feature on the string that produces the cusp, and r is 
the distance between the cusp and the point of observa- 
tion. In natural units Gp/c 2 can be thought of as the 
dimensionless mass per unit length, or tension, of cosmic 
strings. The size L of the feature on the string that pro- 
duces the cusp also determines the low frequency cutoff 
/;. Since L is expected to be cosmological, for example 
the size of a cosmic string loop, the low frequency cutoff 
of detectable radiation is determined by the low frequency 
behavior of the instruments: for the LIGO instruments 
by seismic noise. The high frequency cutoff depends on 
the angle 9 between the line of sight and the direction 
of the cusp. It is given by fh ~ 2c/(8 3 L) and can be 
arbitrarily large (up to the inverse of the light crossing 
time of the width of strings). 

Following [3T] for our templates we take 

r(/)=r 4/3 e(/ h -/)9(/-/0, (2) 

so that h(f) — At(J). We can normalize our tem- 
plates by defining the detector-noise-weighted inner 
product |40j in terms of the two frequency series x(f) 
and y(f) as 

(3) 

Here Sh(f) is the single-sided spectral density defined by 
(n(f)n*(f)) = \8{f - f')S h {f) and n(f) is the Fourier 
transform of the detector noise. We take the inner prod- 
uct of a template with itself to be a 2 = (r|r) and define 
the normalized template f = r/a, so that (f |f ) = 1. 

The calibrated output of an interferometer can be writ- 
ten as 

s(t) = n(t) + h(t), (4) 

where n(t) is the instrumental noise, and h(t) = F + h + (t) 
the gravitational wave signal. The antenna pattern re- 
sponse function to +-polarized gravitational waves, F+, 
is a function of the sky- location of the cusp and the polar- 
ization angle. The signal to noise ratio (SNR) is defined 
in terms of the inner product as p = (s|r). For the case 
of Gaussian noise and in the absence of a signal the SNR 
is a Gaussian variable with zero mean and unit variance. 
In the presence of a signal of amplitude A, the signal to 
noise ratio is a Gaussian random variable, with mean Aa 
and unit variance. Since the average SNR is (p) = Aa, 
for a particular realization of the measured SNR p we 
can identify a signal amplitude 

A = p/a, (5) 

which has an average value (A) = F + B. 

In our search we set // = 50 Hz to be our low frequency 
cutoff. Due to the low frequency behavior of our instru- 
ments, a negligible SNR would be gained by including 
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frequencies lower than 50 Hz. We look for signals with 
high frequency cutoffs fh in the range 75 Hz-2048 Hz. 
The sensitivity of the instruments is such that very lit- 
tle is lost by limiting the search to signals with high fre- 
quency cutoffs above 75 Hz and below 2048 Hz. The only 
template parameter is the high frequency cutoff fh and 
the template bank (the set of templates the determine 
the signals we search for) is constructed iteratively by 
computing the overlap between adjacent templates |5T] . 
The maximum fractional loss of signal to noise is set to 
0.05 and along with the spectrum determines the spac- 
ing between the high frequency cutoffs of the different 
templates. The spectrum Sh(f) is estimated using the 
median-mean method |41j which is fairly robust against 
non-stationarities in the data, including loud simulated 
signal injections. 



C. Trigger generation 

To produce our trigger data we proceed as follows. We 
apply the matched-filter for each template and all possi- 
ble arrival times using fast Fourier transform convolution 
(as described in [31]). This procedure results in a time 
series for the SNR sampled at 4096 Hz for every tem- 
plate. In each of these time series we search for clusters 
of values above the threshold pth — 4 which we identify 
as triggers. For each trigger we determine 

(a) the SNR p of the trigger (the maximum SNR of the 
cluster), 

(b) the peak time of the trigger (the location in time 
of the trigger SNR), 

(c) the start time of the trigger (the first value above 
threshold in the cluster), 

(d) the duration of the trigger (the length of the clus- 
ter), 

(e) the high frequency cutoff of the template, and, 

(f) the amplitude of the trigger A, given by A = p/er, 
where a is the template normalization. 

When several templates result in triggers that occur 
within a time of 0.1 s we select the trigger with the largest 
SNR within that time window. 

We apply this procedure to the data sets of the three 
LIGO interferometers to produce a list of triggers for each 
instrument. 



D. Trigger consistency checks 

To reduce the rate of events unassociated with gravita- 
tional waves (noise induced events) we demand that the 
peak times of triggers in each instrument be coincident 
in time with triggers in the other two instruments. The 



time window used for H1-H2 events is 2ms. This coinci- 
dence window allows for calibration uncertainties as well 
as shifts in the peak times of triggers induced by fluc- 
tuations in the noise. For coincidence between events in 
either of the two Hanford instruments with events in the 
Livingston interferometer a 12ms coincidence window is 
used. This allows for the maximum light travel time be- 
tween sites of 10ms along with calibration uncertainties 
and shifting of the peak location due to noise. We re- 
quire strict triple coincidence: in order for an event in 
one interferometer to survive it must be coincident with 
events in the other two instruments and those two events 
must also be coincident. 

Additionally, we impose a symmetric consistency check 
on the amplitudes of HI and H2 coincident events [5Tj . 
In particular, for HI events we demand that 



\Am-A M2 \ 
A m 



< 



S 

Pm 



(0) 



along with an analogous requirement for H2 events. Here, 
5 is the number of standard deviations of amplitude dif- 
ference we allow and k is an additional fractional dif- 
ference that accounts for other sources of uncertainty 
such as the calibration. We conservatively set S = 3 
and k = 0.5. The purpose of this loose cut on the am- 
plitude of events is to eliminate large SNR events seen 
in one instrument but not in the other. These events 
are not due to gravitational waves but rather to instru- 
mental glitches in the data streams. It is worth pointing 
out that because each instrument has its own antenna 
response factor F + that depends on the orientation and 
direction of the gravitational wave, the H1-H2 amplitude 
consistency test cannot be applied to Hl-Ll or H2-L1. 
Instruments that are not co-located may have different 
values of F + and therefore the measured amplitudes A 
could differ significantly. 

Both time coincidence and amplitude consistency 
checks reduce the rate of events in each instrument from 
about 1 Hz to about 10 //Hz. To simplify the analysis 
in the following we only use the results and statistics for 
the HI events. 



III. FOREGROUND AND BACKGROUND 

To estimate the rate of accidental coincidences that 
survive our consistency checks, the so-called background, 
we time shift the trigger sets relative to one another and 
look for coincident events. We do not perform time shifts 
on the two Hanford trigger sets. Environmental distur- 
bances are known to cause correlations between triggers 
in both interferometers and the effect of time shifting 
the two Hanford trigger sets would be to underestimate 
the background. We therefore take the double coincident 
HI and H2 triggers (which already satisfy the amplitude 
consistency check) and time-shift them relative to the LI 
trigger set. This amounts to treating the two Hanford 
instruments as a single trigger generator, on the same 
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H1 Amplitude (s~ i) Injection Amplitude (s~ I) 



FIG. 1: Plot of the cumulative rate of events as a function 
of the amplitude for both foreground events (filled circles), 
as well as average number of events found in the time shifted 
data (stair steps). The shaded region corresponds to 1-a vari- 
ations measured in the time shifts. 



footing as LI but with a much smaller trigger rate. For 
each trigger in each time shift we then demand the the 
first consistency criterion be satisfied, namely, that each 
Hanford trigger peak be within 12ms of a Livingston trig- 
ger peak. We performed 100 time shifts, with Livingston 
triggers shifted by approximately 1.79 s, the total time 
shift ranging from —89.3 s to 89.3 s. The time shifts 
are sufficiently large that coincident events cannot result 
from gravitational wave bursts from cosmic string cusps. 

Figure [T] summarizes the results of this procedure for 
the HI trigger set. We plot the cumulative rate of events 
for both foreground (unshifted) events (filled circles), as 
well as the average rate of events found in the time shifted 
data (stair-steps) binned in amplitude. The shaded re- 
gion corresponds to 1-a uncertainties computed from the 
variations in the number of events found in the time 
shifted data. 

The loudest HI event has an amplitude of A L = 3.4 x 
10 -20 g -i/3_ There 

are no foreground events which devi- 
ate significantly from the time-slides, and a Kolmogorov- 
Smirnov test confirms that the foreground and back- 
ground distributions are consistent at the 77% confidence 
level. We therefore conclude that no gravitational waves 
have been detected in this search. 



IV. EFFICIENCY 

To determine our sensitivity and construct an upper 
limit we injected over 7400 simulated cusp signals into 
our data set and performed a search identical to the one 
described above. 

The distribution of high frequency cutoffs fh for the 

— 5/3 

injected signals is dN oc f h dfh, appropriate for the 
cusp signals we are seeking (31]. The lowest high fre- 
quency cutoff injected is /* = 75 Hz, coinciding with 
the lowest high frequency cutoff of our templates. The 



FIG. 2: Plot of the detection efficiency e(B) (fraction of events 
detected) as a function of the injected amplitude B. The solid 
curve corresponds to injections with recovered amplitudes A 
larger than our loudest event A L = 3.4 x 10 -20 s -1 / 3 , and 
the dashed curve to all recovered injections. The shaded 
regions indicate our uncertainty in the efficiencies and am- 
plitudes that result from counting, binning and calibration 
systematics. 

amplitudes are distributed logarithmically between B = 
6 x 10~ 21 s -1 / 3 and B = 10~ 17 s -1 / 3 spanning the range 
of detectability. The sources are placed isotropically in 
the sky with sufficient separation in time so as not to un- 
duly bias the spectrum estimate needed to perform the 
matched-filter. 

An injection is found if its peak time lies between the 
start time and the end time of an HI triple coincident 
trigger. We record the recovered amplitude of the injec- 
tion, which is different (typically smaller) than the in- 
jected amplitude because of antenna pattern effects as 
well as noise induced fluctuations. 

The result of our injection run is summarized in Fig. [2] 
We plot the efficiency e(B) (the fraction of injections de- 
tected in the triple coincident HI trigger set) as a func- 
tion of the injected amplitude B for recovered HI am- 
plitudes A greater than our loudest event A L = 3.4 x 
10~ 20 s -1 / 3 (solid line), and for any recovered amplitudes 
(dashed line). The shaded regions indicate our uncertain- 
ties in the efficiencies and amplitudes. The procedure by 
which these curves and their uncertainties are produced 
are described below in more detail. 

A useful measure of our sensitivity is the injected am- 
plitude at which we recover half of our injections. For all 
recovered injections this amplitude is 

B 50% = (5.2±0.9)xl0- 2 V 1 / 3 , (7) 

and for those recovered with amplitudes above our loud- 
est event above it is 

B 5 i o% = (Ll±0.2)xl0- 19 s- 1 / 3 . (8) 

These numbers are consistent with our expectations. 
In [3T] a sensitivity estimate was made for Initial LIGO 
B\lm° ~ 10~ 20 s -1 / 3 . The current search is somewhat 
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less sensitive for two reasons. First, data from the S4 run 
is about a factor of two less sensitive than Initial LIGO. 
Second, we demand coincidence with events in the H2 in- 
terferometer which is about a factor of two less sensitive 
than HI. Together these account for a factor of about 
four leaving only about a 30% discrepancy between the 
rough Initial LIGO sensitivity estimate made in |31] and 
Eq. 0. 

As stated above, the software injections are generated 
randomly, with injected amplitudes that are uniformly 
distributed in their logarithm. Individually, each injec- 
tion is either found or missed. To estimate the probabil- 
ity of injections with a given amplitude being recovered 
we used a sliding window to count the number of software 
injections that were made and recovered within an inter- 
val around that amplitude. The window used was Gaus- 
sian in the logarithm of the injected amplitudes. Choos- 
ing different widths for the window will yield qualitatively 
equivalent but quantitatively different efficiency curves. 
As pointed out above, a useful measure of our sensitivity 
is B 5 oy , and we chose the width of the Gaussian window 
to minimize the uncertainty in B 50 %. 

Three uncertainties are associated with the efficiency 
curve. First, because at each point the value of the ef- 
ficiency has been measured by counting a finite num- 
ber of injections, there is an uncertainty in the efficiency 
attributable to binomial counting fluctuations. Second, 
there is an uncertainty in the amplitude to which a mea- 
surement of the efficiency should be assigned, on account 
of it having been found by counting injections spanning 
a range of amplitudes. Finally, uncertainties in the cal- 
ibration translate into an additional uncertainty in the 
amplitude, on account of the injections from which the 
efficiency was measured having been done at amplitudes 
different from what was intended. The calibration un- 
certainty we use is 11%. This number results from a the 
systematics in the calibration models (5%) and our use 
of time domain calibrated data (10%), which we com- 
bine in quadrature. The counting, amplitude range, and 
calibration uncertainties described above above are com- 
bined in quadrature to produce the shaded regions shown 
in Fig. [2] 



V. PARAMETER SPACE OF COSMIC STRING 
MODELS: CONSTRAINTS AND SENSITIVITY 

For simplicity in this section we will adopt units where 
the speed of light c = 1. The parameter space of cosmic 
string models we need to consider depends on whether 
loops in the cosmic string network are short-lived (life- 
time much smaller than a Hubble time) or long-lived (life- 
time much larger than a Hubble time). This, in turn, 
depends on loop sizes at formation. If their size is given 
by the gravitational back reaction scale [42, 43 , then to 
a good approximation all loops have the same size at for- 
mation and they are short lived. In this case their size 
at formation at cosmic time t can be approximated by 



I = eTG/it, where e < 1 [30] is an unknown parameter 
that depends on the spectrum of perturbations on cos- 
mic strings, and T is a constant related to the lifetime 
of loops, and is measured in simulations to be T ~ 50. 
Recent cosmic string network simulations, however, sug- 
gest loops form at much larger sizes given by the net- 
work dynamics [HI U5]- If this is the case it has been 
shown [B8[ that the regions of parameter space accessible 
to Initial LIGO arc already ruled out by pulsar timing 
experiments. So here we will consider only the first pos- 
sibility, that loop sizes are determined by gravitational 
back reaction and take the size of loops at formation to 
be I = eTGfj,t. 

Unlike field theoretic cosmic strings, cosmic super- 
strings do not always reconnect when they meet. Rather, 
they reconnect with probability p which has a value in the 
range 10 _3 -1 [9]. The reasons for this are (1) that fun- 
damental strings interact probabilistically, and (2) that 
in more than 3 dimensions, even though two strings may 
appear to meet in 3 dimensions they miss each other in 
the extra dimensions. The effect of the decreased recon- 
nection probability is to increase the density of strings by 
a factor inversely proportional to the reconnection prob- 
ability [30J. 

For a point in cosmic string parameter space (G^t, e,p) 
we can use the efficiency curves e(-B) to compute the rate 
of bursts we expect to observe in our instruments, which 
we will refer to as the effective rate 7. It is given by the 
integral [3T] 



' e{B) dR(B;G^e,p) d ^ 
dB 



where B ~ Gjjd 2 / 3 /r is the optimally oriented amplitude 
(i.e. the amplitude of events excluding antenna pattern 
effects), e(B) is the efficiency of detecting events at an 
amplitude B and dR(B; G/i, e,p) is the cosmological rate 
of events with amplitude in the interval B and B + dB. 
We have take the size of the feature that produces the 
cusp to be the size of the loop I. 

Since we are considering loops that are small when 
they are formed they are also short lived, and at a given 
rcdshift they are all of essentially the same size. As a 
result the amplitude of burst events from a given rcdshift 
is the same. In this case, rather than Eq. ^ it is easier 
to evaluate, 



7(G/i,£,p) 



€{z) dR(z;G^e,p) dz _ (1Q) 



Here dR(z) the rate of bursts originating at redshifts in 
the interval between z and z + dz. The rate is given by 
Eq. (59) of [SI] 

dR „ N c {g 2 UH^)-^ _ 14/3 5/3 
Tz = H ° 2a^pTGv + 



x 6(1 - 6 m (z,f*,aH VtO)), 



(11) 



where Hq is the present value of the Hubble parameter, 
iV c is the average number of cusps per loop oscillation; 
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g2 is an ignorance constant that absorbs the unknown 
fraction of the loop length that contributes to the cusp 
and other factors of O(l); /* is the lowest high frequency 
cutoff of the bursts we are interested in detecting; a — 
eTG/i is the loop formation size in units of the cosmic 
time; and 9 m = [52(1 + z)/*^] ^ 3 is the maximum angle 
a cusp and the line of sight can subtend and still produce 
a burst with high frequency cutoff /*. The 0- function 
removes events that don't have the form of Eq. 0. Two 
dimensionless cosmological functions enter the expression 
for the rate of events: ift(z) which relates the cosmic time 
t and the redshift via t = H^ip^z), and fv( z ) which 
determines the proper volume element at a redshift z 
through dV(z) = H ~ 3 tpv(z)dz (see Appendix A of [3T]). 
For details on the derivation of this expression see [2"51 

nana no- 

The efficiency e(z) is the fraction of events we detect 
from a redshift z. We compute this quantity starting 
from our measured efficiency as a function of the ampli- 
tude B, using Eq. (60) of [51] 

£'\z) = BH~ 1/3 

where g\ is an ignorance constant that absorbs the un- 
known fraction of the loop length that contributes to 
the cusp and factors of 0(1) (different from those of c/2, 
see [31| for details), and ip r (z) is a dimensionless cosmo- 
logical function that relates the proper di stan ce r to the 
redshift via r — 

ff ~~VO)- Solving Eq. Q for z gives 
the redshift from which a burst of amplitude B originates. 
Thus for each amplitude in our efficiency curve we can 
determine the corresponding redshift and construct e(z). 

We have computed the cosmological functions ip r (z), 
ift(z), and (fiv(z), as described in Appendix A of [51] . 
using the latest consensus cosmological parameters mea- 
sured by WMAP [46]. Specifically, we used a present 
day Hubble parameter of Hq = 70.1km s _1 Mpc _1 
'2. '27 / 10~ 18 s _1 , and densities relative to the critical den- 
sity f2 m = 0.279 for matter, $2l r = 8.5 x 10~ 5 for radia- 
tion, and J1a = 1 — — ^l r — 0.721 for the cosmological 
constant. 

We scan the parameter space of the small loop models 
by varying the reconnection probability p, the dimension- 
less string tension G/x, and the size of loops parametrized 
by e. To construct a 90% upper limit we use the loudest 
event statistic [57]. For each point in cosmic string pa- 
rameter space we compute the effective rate 7 using the 
solid black curve in Fig. [2] (recovered injections with am- 
plitudes above our loudest event). We then compare this 
rate with 790% = 2.303/T, the rate for which 90% of the 
time we would have seen at least one event in a Poisson 
process if we observed for a time T. If for a point in pa- 
rameter space the effective rate 7(G/i, e,p) exceeds 790% 
then we say those parameters are ruled out at the 90% 
level, in the sense that cosmic string models with those 
parameters would have produced an event with ampli- 
tude larger than our loudest event 90% of the time. We 
can also estimate our sensitivity by computing 7 using 



the dashed curve in Fig. [2] (recovered injections with any 
amplitude) and then compare it to \/T. The latter tells 
us what models would have resulted in at least one event 
(on average) surviving our pipeline, a useful measure of 
our sensitivity. For our search in triple coincident S4 data 
T = 14.9 days. 

Figure [3] shows the results of our analysis. We have set 
N c = gx = gi = 1 for convenience. Each of the four pan- 
els corresponds to a value of the reconnection probability 
p with the loop size parameter s on the y-axis, and the di- 
mensionless string tension Gfj, on the cc-axis. Regions to 
the right of the red curves are excluded at the 90% level 
by our analysis. The dotted curves show the uncertainty 
in the areas of parameter space excluded that result from 
uncertainties in the efficiency The light and dark gray 
shaded areas in each of the four panels are regions of pa- 
rameter space very unlikely to result in a cosmic string 
cusp event detected in S4: a cosmic string network with 
model parameters in these regions would result in less 
than one event (on average) surviving our pipeline. The 
lower boundary of these areas was computed using the 
efficiency for all injections, the dashed curve shown in 
Fig.[2j and comparing the effective rate to 1 /T. The dark 
gray regions show regions of cosmic string model param- 
eter space unlikely to result in a cosmic string cusp being 
detected in a year long search of Initial LIGO data: on 
average such models would result in fewer than one event 
being detected. The lower boundary curve of the dark 
gray regions was computed with the Initial LIGO sensi- 
tivity estimate B^° w 10~ 20 s~ 1/3 [5T] and assuming 
a year of observation time. 

Due to the sensitivity and duration of the S4 run the 
90% limits we have placed on the parameter space of 
cosmic string models are not as constraining as the indi- 
rect bounds due to Big Bang Nucleosynthesis [55]. An- 
other current gravitational wave bound comes from pul- 
sar timing observations. Due to their sensitivity at very 
low frequencies pulsar timing bounds for loop sizes given 
by gravitational back reaction constrain an independent 
portion of the cosmic string parameter space [38] . Future 
analysis of data from the fifth science run, however, a fac- 
tor of two more sensitive and with a year of triple coin- 
cident data, will be sufficiently sensitive to surpass these 
limits in large areas of cosmic string parameter space. 



VI. SUMMARY 

We have performed a search for bursts from cosmic 
string cusps in 14.9 days of triple coincident data from 
LIGO's fourth science run. The gravitational waveforms 
of cosmic string cusps are known and matched-filters pro- 
vide the optimal means of extracting such signals from 
noisy data. We constructed a template bank and gener- 
ated a set of triggers for each of the three LIGO interfer- 
ometers. To reduce the rate of accidentals we demanded 
the resulting triggers satisfy two consistency criteria: (i) 
time coincidence between events in the three instruments, 



9 



1 

10- 2 

10-" 
os io -6 

10" 8 
10 -10 



p=10 



p=10 




1 

IO" 2 

IO" 4 
a 1 g-B 

10" 8 
10 -10 

10" 12 



10"I2 I i i i i S i ' i "-V- 1 

10" 12 10" 11 10" 10 10" 9 10" 8 10" 7 10" 6 10" 5 10" 4 



p=10 



















I 




\ 


\ 






10" 12 10" 11 10" 10 10" 9 10" 8 10" 7 10" 6 10" 5 10" 4 

p=1 



10" ; 



10"' 



10"' 



10" 1 



10" 1 



10" 12 10" 11 10" 10 10" 9 10" 8 10" 7 10" 6 10" 5 10" 4 
Gn 



10" 



10" 12 10" 11 10" 10 10" 9 10" 8 10" 7 10" 6 10" 5 10" 4 



FIG. 3: Plot of the sensitivity and upper limit results of our analysis. Each of the four panels corresponds to a value of the 
reconnection probability p with j/-axes e and :r-axes Gp. Areas to the right of the red curves show the regions excluded at 
the 90% level by our analysis. The dotted curves indicate the uncertainty in the areas of parameter space excluded that arise 
from uncertainties in the efficiency. The light and dark gray shaded areas are regions of parameter space unlikely to result 
in a cosmic string cusp event detected in S4: a cosmic string network with model parameters in these regions would result in 
less than one event (on average) surviving our pipeline. The lower boundary of these shaded areas was computed using the 
efficiency for all recovered injections, the dashed curve shown in Fig. [2] The dark gray regions show regions of parameter space 
unlikely to result in a cosmic string cusp being detected in a year long search of Initial LIGO data. The lower boundary of the 
dark gray areas was computed with the Initial LIGO sensitivity estimate -650%° ~ 10 -20 s -1 ^ 3 |31j . 



and (ii) amplitude consistency between events in the two 
Hanford interferometers HI and H2. The latter check is 
possible because the Hanford instruments are co-aligned. 
The effect of the consistency criteria is to reduce the trig- 
ger rate in each of the instruments from about 1 Hz to 
about 10 pHz. To estimate our background, the rate of 
accidentals, we performed 100 time-slides on the data, 
with a time step much larger than the duration of our 
signals. Comparing our background estimate with our 
foreground, we conclude no gravitational waves have been 
found in this search. 

To estimate the sensitivity of the search and place con- 
straints on the parameter space of cosmic strings we per- 
formed several thousand simulated signal injections into 



our data streams and attempted to recover them. The 
simulated signal parameters used are consistent with our 
expectations for the cosmic string population. The in- 
jections were used to compute the efficiency, the fraction 
of events that we detect in a range of amplitudes. The 
efficiency curves can be convolved with the cosmological 
rate of events to compute the so-called effective rate: the 
rate of detectable events in our search. The effective rate 
folds together the properties of the population, such as 
the distribution of sources in the sky, with the sensitivity 
of the detectors and the analysis pipeline. We found our 
estimate for the sensitivity of this search to be consistent 
with our expectations, given that we used S4 data and 
included H2 in the analysis. Using the loudest event in 
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our foreground, we placed a frequentist 90% upper limit 
on the effective rate, which we in turn used to constrain 
the parameter space of cosmic strings models. Unfortu- 
nately, the sensitivity of this search does not allow us 
to place constraints as tight as the indirect bounds from 
Big Bang Nucleosynthesis. However, analyses using data 
from future LIGO runs is expected to surpass these limits 
for large areas of cosmic string parameter space. 
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